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TITLE 

TYPE 
DURATION 



ACCURACY 
DESCRIPTION 



DATA TAPE PREPARATION 



Complete Linear fcfetrix Equation Solver and General 

Matr ix Inversion 

Complete program 

n = order of matrix A, m = number of columns of matrix B. 

a) reading in of matrix A and Bs read- in speed 
(250 characters ner second) 

b) computing time} a /5OO (m = l) 

1 -11 + kn. « ,3n + .Oln (n = m) 

c) punching of matrix Xj punch speed (60 characters 
per second) 

A function of the order and conditioning of the matrices . 
This routine solves the linear: matrix equation AX = B, 
for matrices A and B satisfying the conditions 1 

a) A is non-singular and of size (nxn) 

h) B is of size (nxm) 
-where the magnitude of n and m is governed "by the 
inequality s 

nm + n/2(n + I) + (n + m) < 842 
The solution X, of size (nxm) > is punched out by 
successive columns * In the special case when B = I 
(the identity matrix) - t we obtain the inverse since 
the solution is X - A* . Here n = m and. the limits 
on the size of A are n < 22 * 

A distinction is to be made between the cases when 
B = I (i»e* when we wish to invert A) and B £ I 
(i.e. when we are in essence solving m sets of n 
simultaneous equations in n uiLknowns ) « 

1) B = I Punch the matrix A, row by row, adhering 
to the following % 

a ) Each row is ended by J 4 

b) Follow the final J with a sexadecimal 
character 1, 2, ..., K, S indicating 
the number of digits desired in the 
punching of X. 



■2* 



2) B £ It Bsach the augmented matrix [A,B] row by 
row, such that i 
a ) Each of the respective rows of A and B 

are ended by N» 
h) Follow" the final IT of the last row of 
B "by a sexadecimal character 1, 2, . .,, 
K, S indicating the number of digits 
desired in the punching of the output X. 
j5) In either eases 

a ) All rows, must he so scaled that no 
element is > l/2# 

b) Ikeh element must he preceded by a 
plus or minus sign and may contain up 
t©' 12 figures »■ (indicate a plus sign 



for zero)* 


X 


2 





A typical input tape for A =* 


3 


«4 1 




-2 


5 6 


. 


r l 


2 




and B = 





-1 




would appear as follows g 


^2 







+1 +2 + I +1 +2 


IT 




+3 -4 +1 B" + «1 


H 




-02 +05 +06 H" -02-8- 


ST K 







and will print out to 10 places rounded off. 
If B were I we would punch instead? 

+01 +02 + J 

+03 -0+ +01 J 

*02 +05 +06 J K 
Notes Concern must be given to the scaling possibilities 
offered by each case* In the former case* each row can 
be considered independently of the others 3 remembering, 
however, that corresponding rows of A and B must be 
treated as a unit* In the latter case, scaling one 
row forces one to scale each row. 



METHOD OF USE 



DATA. OUTPUT 



MATHEMATICAL METHOD 



1 ) Bead in the routine until it stops « 

2) Flace the data tape into the reader and start "by 
raising the black switch* 

5) The computation is then carried out and subsequently 
X will "be punched out* The machine will then stop on 
a transfer order (2^- 52) and is then ready to accept 
another equation for solution upon raising the 
"black switch* !o concern need "be given to successive 
values of matrix sizes - 
Each column of X will he punched out in succession* 
followed by a scaling factor and the character I* This 
scaling factor may he different for each column* It 
indicates the position of the decimal point as lying 
after the eolumn in which the factor 1 appears* A 
typical output for a 3 x 3 matrix may look like the 



following? 








+1000 








-0500 
+2500 
+1000N 
-2000 




1 


~2Q 


+0050 
-2500 
+0100N 


representing 


-5 
2.5 


.-5 
-25 


+0500 
+0100 








+0050 
+1000H 









.3 
.1 
.05 



Successive calculations are separated by three carriage 
returns and line feeds. 

Motes In case we are inverting A(i*e.j B = i), then 

n+1 
the output matrix must he multiplied hy 10 if the 

original matrix A was scaled hy 10 * 

1) The matrix A is upper triangularized by a series 
of elementary row operations j, TTE. * during which 
time the same operations are performed upon the 
augmented matrix B* obtaining thereby (TTE. ) B = B 1 ' «■ 

2) Successive columns of B T are then considered and 
by the back substitution method we solve m sets of 
n simultaneous equations in n unknowns* The m 



A. 



column vectors comprising this solution form 
the matrix X* 

5) For the case B - I, it follows that X = A ' . 

k) In the event a zero appears on the diagonal of 
the upper triangular matrix, the machine detects 
the singularity of A and prints out F, indicating 
that there is no solution. It is then ready to 
proceed, to the next problem* 

5) Appearance of a zero on the diagonal of the upper 
triangularized matrix is a necessary condition for 
singularity of A* This is not true , however, for 
the machine representation, due to accumulating 
round off errors. In the process of upper triangular- 
izing A* automatic scaling is performed so that 

the absolute value of the largest element in each 
row of the augmented matrix lies in the interval 
(l/2, l/43. This destroys the possibility of 
comparing the relative size of the diagonal 
elements upon the completion of "the calculation, 
in order to obtain an indication of the rank of 
the matrix # It is, therefore, impossible to 
guarantee detection of a nearly singular matrix. 
However, it should "be noted that the size of the 
sealing factors appearing for each column of X 
will indicate the conditioning of A* In the event 
that scaling gets beyond the capacity of the 
registers of Illiac the calculation is stopped 
and F is printed out • The machine is then ready for 
another calculation* 

6) In addition to the precautions of this routine, it 
is advisable, in questionable instances, to form 
the product of A and the X f computed, by this routine^ 
and note the approximation to the matrix B* 
Another alternative would be to employ Routine M 12 
by which one can find the determinant of A and can 
inspect the elements on the diagonal of the upper 



HOTES 



triangularizied A and obtain the rank, 

1) The triaftgtilarizatiQn process and back substitution 
is a modification of Koutine L 1. 

2) Since (A )* = (A ) e© e eaa either punch rows 
or columns depending on whether the input 



of A 



>T 



matrix is A orA respectively* 
3 ) For the case B = I, the identity matrix is 

automatically augmented within the routine. 
k) A sample inversion of a matrix of order 22 



required approximately 35 seconds computing time 

-11 
Errors were less than 10 * 

5 ) A sample solution of a system of 38 equations 

in 38 unknowns with random integer coefficients 

required approximately 100 seconds computing 

-10 
time* Errors were less than 10 . 
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CODED BY 



AHPR0YED BT 




loc&tcou 



6 



7 



10 



11 



12 



13 



l* 



15 



NOTES 
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Routine X 1 

00 8K 
Routine N 2 

00 3^K 
Routine F 2 

00 52K . 

4l ?r 

41 6? 
50 183F' 
50 1L 
26 Sf 
ito F 

L5 16F 
J+6 17L 
46 83L 
46 91L 
10 20F 

42 82L 
42 9% 

. 42 10L 
42 126L 
4-2 12% 
LO 46L 

42 4F 
42 3F 
L3 F 
30 17L 

41 (y)F 
F5 10L 

42 ItL 
LO 3F 
LO 126L 

. 32 10L 
L5 126L 
L4 TF 
42 15L 
L5 125L 
4o (y+r)F 



6* 11* 



14' 



Decimal Order Input 

Input a Sequence of -Decimal Fractions 

Print (A) with or without Sign to n Places 

Clear row counter 
Clear counter 

— Read in row r of A 

Store or 1 

183 + n = y (left address ) 



Set address of y 



y - 183 
Store m = n 
Store n 

Test for inversion or solving AX = B 



— Clear y to (y + n - l) 



y + r 
1/10 



Augment unit matrix 



LOCATION 


OHEER | NOTES MGE 2 | 


16 


22 20L 




i 




00 t 




1 

Waste 


17 


50 (y)F 

50 1?L 


3' 




— Eea& ±m row r of B if : B ^ T 


18 


26 ar 










L5 l6F 




y + m - t 


19 


10 20F 








LO 126L 




t « y = m 


20 


42 4f 
U~> 46% 




Store m 

j 


21 


42 26L 




i 
j 




L5 4F 






\ 


22 


I»4 126L 

42 48L 






1 

i 

j 


23 


42 2JL 
42 1281, 






rt 1 
— I03 + n + m = t j 


24 


00 20F 
46 127L 






| 


25 


22 25L 
L3 TF 






Waste | 

| 


26 


32 27L 






Cfemse Itb. row to inter cbaoge -with 




4l ( )F 


21 




Tirttml it& row 


27 


26 28L 










41 (t)F 


23 






28 


L5 48l 
4-2 38L 








29 


k2 4lL 
42 44L 








30 


42 ^4L 
L5 46l 








31 


l4 6f 
42 4-7L 




183 + i = x 




32 


42 51L 
42 68L 




* 


— Set addresses 


33 


42 TOL 
42 42L 






1 



LOCATION 


ORDER 


NOTES PAGE 3 


34 


00 20F 
46 4lL 








35 


46 39L 
46 45L 








36 


46 51L 










46 53L 






37 


L5 123L 
46 46% 






— Prepare for interchange of rows 


38 


46 49L 








39 


L3 (t)F 
L6 (x)f 
32 42L 


28 « 

35" 


Test size of leading elements 


40 


47 48L 




No row interchange 




50 7F 




Approximate zero in Q 


4i 


L5 (x)F 


34 J 






66 (t)F 


29 




42 


26 46L 








50 (x)F 


33 ' 




43 


S3 F 
32 45L 




- 


44 


50 124L 




1 .. 2-39 




75 (t)F 


29 » 




45 


66 (x)F 


35' 






47 4% 




Row interenange 


46 


4l 5F 




Clear counter 




SI 183F 




Address a constant 


47 


40 2F 








L5 (x)F 


31 l 58 1 






48 


40 ( )F 
L5 (t)F 


37* 4o 
22' 56 






^9 

1 
t 
i 

t 


40 ( )F 

50 2F 


38 45 * 






i 50 

t 


7J IF 
L4 F 








I 51 


40 (x)F 


36 57 T 






1 


L3 (x)F 


32 57* 


~I 





M 13 



LOCATION 


ORDER 


_ _ — — 1 

NOTES PAGE 4 M 13 


52 


L6 5F 


|H(5)| - |N(x)| 


Element of 


Linearly 




36 5^L 






largest ah-* 


combine i 

1 


53 


L7 (x)F 


36' 58 




solute value 


successive! 




. 40 5F 






in row i 


rows so j 

j 


5^ 


L5 IF 






as to get 




40 (t)F 


30 55 1 




zeros ' 


55 


F5 5^L 
42 jkL 






"below the j 
diagonal j 


56 


42 48L 
L5 51L 






of A 

i 


57 


L4 123L 
40 51L 






! 
i 


! 58 

i 

1 
j 


46 53L 
42 47L 








1 59 

1 


LO 127L 
32 kjL 








i 6o 


L5 74L 
40 69L 




Reset order 


1 61 


L3 5F 
36 75L 




If zero do mot rescale 


j 62 

t . 

63 


li 5F 
32 65L 
L5 70L 
46 69L 




If absolute -value of any one element 
of row is > l/2 prepare order to 
scale down 


64 


26 68l 






: 


F5 69L 






65 


42 69L 




If absolute value of largest element 


| 66 

1 

t 

i 67 


F5 69L 
42 69L 
L5 5F 
00 IF 




< l/s prepare to scale -up rows so 
that this element is in the interval 
(1/2, l/k] 


i 


4o 5F 






; 68 


50 7F 
L5 (x)F 


72 32 r 


Put ~ zero in Q 




i 69 

f 

j 1 


10 (1)F 
00 (l)F 


63* 60 * 
66 65 60 « 


- 





LOCATION 



70 

71 
72 

73 
7^ 
75 
76 
77 
78 

79 
80 
81 
82 
83 

eh 
85 

86 

87 
88 



ORDER 



50 2F 
40 (x)F 
F5 70L 
42 7QL 
42 68L 
LO 128L 
36 68L 
26 75L 
10 IF 
00 IF 
F5 6F 

4o 6f 

L5 7F 
LO 6F 
36 28L 

Fp 7^ 
40 7F 
LO 3F 
36 80L 
22 L 
in 5F 
81 4f 
00 20F 
46 115L 
L5 125L 

40 (y)F 
L3 (y)F 
32 121L 

^1 6f 

41 7F 
L5 ^6L 

42 114L 
L5 48L 
42 96L 
LO 4F 
L4 5F 
42 91L 
L5 96L 



NOTES 
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71* 33 



Waste s address a constant 



5' 

4 



Test 



Rescale row i 



Constant 



Determine if row i must haye further 
eliminations (done i times) 



— Count number of rows 



Repeat for next row 
Clear counter 



Read in character to determine number 
of digits to punch 



Set scaling factor 



Terminate calculation if scaling < 2 



Clear counters 



-39 



Set addresses 



LOCAITIOI 


ORDER 


TOES 


PAGE 6 


89 


LO 4F 
42 96L 








90 


L5 91L 
46 96L 








91 


50 (y)F 
71 ( )F 


4' 

88 94 








92 


40 F 
L5 91L 










93 


FO 7F 
LO 4f 










94 


42 91L 
22 100L 






— Calculate 


n 

Z a..x ± 

4-1+1 1J * 


95 


22 95L 
S5 F 




Waste 






96 


50 (y)F 
74 ( )F 


90* 101 * 

86 l 89* 101* 








97 


L4 F 
40 F 










98 


LL F 
32 100L 




> 1/2? 




99 


50 125L 










7J (y)F 


6 


-Hescale and start again 




100 


L5 9& 








101 


LO 123L 








102 


4o 96L 

42 106L 
42 109L 




"Set addresses 




103 


46 TUT. 

F5 6F 








104 


40 6f 
L5 W 




- Count (11 - i + l) times 


for roTsr 1 


105 

i 
» 


LO 6F 
32 95L 








I 106 

1 
1 
l 

i 


41 6y 
L3 (a j[i )F 


102 


Reset 





Loc&noff 


ORDER 








& PAGE 7 




mm 


107 


32 1211, 


1 End if zero on diagonal 




L6 F 




Test if division is proper 


108 


36 9% 








26 109L 




Waste 


109 


id F 








66 (a..)F 

v 11 


102' 




110 


22 110L 
SI F 




Waste 


111 


i^O (y-d)F 


103 






F5 7F 








112 


ho 7F 
LO 3F 






— Corart n rows 


113 


36 ll^L 










22 88L 




Repeat 


nij- 


92 131F 
L5 (x)F 


85 * 117 


Line feed 




115 


50 ( )F 


81* 










50 11% 






- To eater P 2 


Ptmcn out eoliama j 


116 


26 3^F 
F5 ll^L 




"■ 


of X 


117 


k-2 ±Xkh 
FO 12% 








118 


36 ll^L 




— 






92 770F 




Bmen U 


119 


F5 5F 
kO 5F 








120 


LO k¥ 
36 122L 






— Count m coluions 


121 


2.6 82L 




Repeat 




92 898F 




Funcn F 


j 122 


92 139F 




!Ehree line feeds "between problems 


1 


2k L 




Ready for next calculation 


! 123 

! 

i 
1 
1 


00 IF 
00 IF 






! 12^ 


7L 4095F 
LL 4095F 




1 - 2" 39 



LOCATION 
125 

126 

127 

128 

129 



ORDER 

00 F 

00 1000 0000 0000 J 

B6 17L 



NOTES 
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hi (y)F 
NO (t)F 
L3 F 
JO 2F 
kO (t)F 
12 131F 
L5 (y)F 

2il- 52N 



1/10 



T 

2 V 



23' 



V 



End constants 



